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Abstract. While kinematical modelling of the low PN velocity dispersions observed in the outer regions of elliptical galaxies 
suggest a lack of dark matter around these galaxies, we report on an analysis of a suite of A-body simulations (with gas) of 
major mergers of spiral galaxies embedded in dark matter halos, and find that the outer velocity dispersions are as low as 
observed for the PNe. The inconsistency between our dynamical modelling and previous kinematical modelling is caused by 
very radial stellar orbits and projection effects when viewing face-on oblate ellipticals. Our simulations (weakly) suggest the 
youth of PNe around ellipticals, and we propose that the universality of the PN luminosity function may be explained if the 
bright PNe in ellipticals are formed after the regular accretion of very low mass gas-rich galaxies. 


INTRODUCTION 

Although it is generally accepted that spiral galaxies are 
embedded in dark matter halos, the presence and amount 
of dark matter in elliptical galaxies is still uncertain. 

The classical way to infer the mass distribution of a 
gravitational system is by analyzing its internal kine¬ 
matics. For an unordered system such as an elliptical 
galaxy, the more mass there is, the higher the velocity 
dispersions. Mathematically, this is expressed through 
the Jeans equation VP = — v VO, where P, O and v are 
the dynamical pressure, gravitational potential and (num¬ 
ber or luminosity) density of the tracer, which in spheri¬ 
cal symmetry becomes: 

= ( 1 ) 

dr r r^ 

where M is the mass, a> the radial velocity dispersion, 
j3 = 1 — of/(2 (7^) the velocity anisotropy (zero for 
isotropic systems) and o, the tangential velocity disper¬ 
sion. The modelling is complicated by projection effects 
and by the intrinsic mass/anisotropy degeneracy in the 
Jeans equation. 

Alas, stellar velocity dispersions in elliptical galaxies 
are difficult to measure (through absorption-line spec¬ 
troscopy) beyond 2 half-light (effective) radii (2Re). 
Planetary nebulae (hereafter, PNe) happen to provide a 
discrete indicator of the gravitational potential that ex¬ 
tends up to 5 Re, and therefore helps constrain the mass 
prohles of ellipticals. Measured velocity dispersions of 
PNe around ellipticals turn out to be low and decreas¬ 
ing with radius (Ciardullo et al. 1993; Mendez et al. 


2001; Romanowsky et al. 2003, but see Teodorescu, in 
these proceedings), and rehned kinematical analyses lead 
to a lack of dark matter around elliptical galaxies (Ro¬ 
manowsky et al.). This conclusion is at odds with the 
predictions of the standard hierarchical scenario in which 
elliptical galaxies are the products of mergers of spi¬ 
ral galaxies of comparable mass (Springel et al. 2001). 
Do the low PN velocity dispersions force us to revise 
our general understanding of the formation of elliptical 
galaxies? In this contribution, we gain considerable in¬ 
sight through an analysis of dynamical A-body simula¬ 
tions of galaxy mergers. 


THE LINE-OF-SIGHT VELOCITY 
DISPERSION DISPERSION PROFILES 
OF MERGER REMNANTS 

We use a series of 10 simulations (Cox et al. 2005) 
of equal mass mergers of disk-rbulgeH-gasH-dark matter 
halo spiral galaxies, with different orbital parameters, 
bulge/disk ratio, gas fraction and feedback efficiency. 

Figure 1 shows snapshots of one of the mergers, illus¬ 
trating the tidal tails and bridge formed during the closest 
approach, after which the galaxies separate and merge at 
second passage, while the tidal tails expand outwards. 
The hnal remnant, viewed here edge-on is oblate. 

Figure 2 shows the line-of-sight velocity dispersion of 
the 10 merger remnants, together with the stellar and PNe 
dispersions of Romanowsky et al. and Mendez et al.. In¬ 
terestingly, the stellar line-of-sight velocity dispersions 
of the dark matter-embedded merger remnants are as low 
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FIGURE 1. Snapshots of the merger of equal mass spiral 
galaxies (see Cox et al. 2005). The lighter zones in the central 
regions of the galaxies represents the young stars formed during 
the simulation, while the remaining smooth greyscale are the 
old stars. 



FIGURE 2. Line-of-sight velocity dispersion vs. projected 
radius, normalized to the effective radius (from Dekel et al. 
2005). The merger remnants are viewed from three orthogo¬ 
nal directions at two late times after the merger. Lower solid, 
dashed, and upper solid curves show the dispersions from all 
stars, young (< 3 Gyr) stars, and dark matter, respectively. The 
IcT scatter is marked by shaded or hashed areas or a thick 
bar. For comparison, are shown the stellar {crosses, Jedrze- 
jewski and Schechter 1989; Binney et al. 1990; Bender et al. 
1994; Statler and Smecker-Hane 1999) and PN (circles, Men¬ 
dez et al.; Romanowsky et al.) velocity dispersions of galaxies 
NGC 821, 3379, 4494 and 4697. The errors are 1 a. Thin curves 
refer to the earlier models (Romanowsky et al.) with (upper) 
and without (lower) dark matter. 


as observed, and surprisingly, consistent with the predic¬ 
tions of no dark matter of Romanowsky et al.. Also, the 
surface density prohles of the simulated merger remnants 
match almost perfectly the observed surface brightness 
profiles (Dekel et al. 2005). In other words, the simulated 



r/Re„ 


FIGURE 3. Radial (3D) profiles (Dekel et al. 2005) of 
velocity anisotropy, for the dark matter (nearly constant solid), 
young stars (dashed), old stars (dotted) and all stars (rising 
solid curves), respectively. 


merger remnants match the observations, while they are 
formed through the mergers of spiral galaxies embedded 
in dark matter halos, which is in contradiction with the 
results from kinematical modelling (Romanowsky et al.) 
that point to a lack of dark matter in elliptical galaxies. 


WHY DID KINEMATICAL AND 
DYNAMICAL MODELLING PRODUCE 
INCONSISTENT RESULTS? 

One can criticize the relevance of these simulations to 
understanding the kinematics of PNe around elliptical 
galaxies as these are not thought to be subject to recent 
(few Gyr) mergers of equal mass spiral galaxies. How¬ 
ever, the low velocity dispersions are also found in merg¬ 
ers of spirals with mass ratio 3 (which are much more fre¬ 
quent) and these low velocity dispersions are long lived 
if one lets the remnant evolve for many Gyr more. 

If the merger simulations analyzed by Dekel et al. 
are relevant, then how can dark matter-rich spirals make 
elliptical-like merger remnants similar to the observed 
elliptical galaxies, while kinematical modelling of the 
PN motions around ellipticals imply little or no dark 
matter (Romanowsky et al. 2003)? 

There are several reasons for the inconsistency be¬ 
tween the kinematical and dynamical analyses. First, the 
orbits of the outer stars turn out to be quite radial, with 
P « 0.5 (Fig. 3 from Dekel et al. and also Dommguez- 
Tenreiro et al. 2004; Abadi et al. 2005), which decreases 
the observed velocity dispersions, since the dispersion 
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FIGURE 4. Mass-to-light ratio at the virial radius, relative 
to that derived from isotropic NFW modelling, as a function 
of observed line-of-sight stellar velocity dispersion at 2 {right) 
and 5 (left)Re. From Mamon and Lokas (2005). 


FIGURE 5. Ratio of Jeans mass (i.e., the mass from the sta¬ 
tionary spherical Jeans equation [1] — with the velocity disper¬ 
sions replaced by rms velocities to take streaming motions into 
account), to true mass averaged over the 10 simulated merger 
remnants, and in spherical shells. 


at a fixed projected radius depends (Mamon and Lokas 
2005) on the mass profile outside that radius and the 
overall anisotropy profile. This radial stellar anisotropy 
arises because the stars selected to lie far from the cen¬ 
ter in the remnant must have come from the inner re¬ 
gions right after the closest approach, and were thrown 
out on elongated orbits in tidal tails (see Fig. 1). This 
is illustrated in Fig. 4, which shows that for the stel¬ 
lar anisotropy of the merger simulations (upper curves: 
‘Nav04-0.18’) the mass inferred at the virial radius with 
the recent representation (Navarro et al. 2004) of ACDM 
density profiles is 2.4 times what would be inferred for 
an isotropic “standard” NFW (Navarro et al. 1996) dark 
matter model (Mamon and Lokas 2005). Although the 
best fit models found in the kinematical analysis are 
moderately radially anisotropic (Romanowsky 2005, pri¬ 
vate communication), our equal mass merger leads to 
even stronger radial anisotropy, and unequal mass major 
mergers lead to even stronger anisotropy. 

Second, the galaxy modelled in detail by Ro¬ 
manowsky et al., NGC 3379, appears circular and 
its PN system shows no rotation (Ciardullo et al.; Ro¬ 
manowsky et al.) and may well be an oblate elliptical 
viewed face-on (Statler and Smecker-Hane 1999), which 
would lead to roughly 10% lower velocity dispersions 
(Dekel et al.). 

Moreover, the galaxy may not be stationary, contrary 
to the assumptions of all kinematical modelling. How¬ 
ever, Fig. 5 indicates that the spherical stationary Jeans 


equation (1) recovers well the true mass profile, so that 
the merger remnants can be modelled as stationary and 
spherical systems (at the meeting, we incorrectly re¬ 
ported a 40% deficit in the mass estimated by the sta¬ 
tionary spherical Jeans equation). 


ARE THE PNE OBSERVED IN 
ELLIPTICAL GALAXIES YOUNG? 

There are several reasons to believe that the PNe ob¬ 
served in elliptical galaxies are young objects, although 
each explanation comes with its caveat: 

1. in the latest model of PN evolution (Marigo et al. 
2004), the brightest PNe are 1 or so Gyr old, com¬ 
parable to our young stars — however, the univer¬ 
sality of the bright-end of the PN luminosity func¬ 
tion (Ciardullo et al. 1989), between old ellipticals 
and young spirals, suggests otherwise (see below); 

2. the young stars in the simulations formed during the 
merger from the gas particles, have lower velocity 
dispersions than the older stars (dashed vs. solid 
curves in Fig. 2) — however, the effect is of order of 
10%, and the old stars fit the observed PNe velocity 
dispersions as well as the young stars; 

3. the reduced kurtoses, hn, of the line-of-sight veloc¬ 
ity distribution of the young stars are larger than for 
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FIGURE 6. Outer vs. inner Gauss-Hermite h/i coefficient 
for simulated young stars and for PNe in observed galaxies 
(the ellipses represent confidence intervals obtained using boot¬ 
straps). The shaded region indicates the h 4 values obtained 
from spectral absorption lines of elliptical galaxies. 


the old stars and more similar to the /14 parameters 
measured in the two galaxies with publicly available 
PN velocity data, while the h 4 values for the old 
stars are closer to the measurements from spectral 
absorption lines (see Fig. 6 ) — however, the match 
is not excellent. 

None of these reasons is compelling. One clever rec¬ 
onciliation (Ciardullo et al. 2005 and Ciardullo, in these 
proceedings) between the observed universality of the 
PNLF between young spirals and old elliptical galaxies 
and the recent Marigo et al. model linking the bright- 
end PNLF with the age of the PNe is to suppose that 
the progenitors of the bright PNe in ellipticals are binary 
systems (blue stragglers) of fairly old stellar mass stars 
rather than much younger 2Mq single stars. 

Alternatively, it may be possible to reconciliate PNLF 
universality and the Marigo et al. PN model if one sup¬ 
poses that the bright PNe in ellipticals are caused by re¬ 
cent (1 Gyr) accretion of very low mass galaxies that go 
through a star burst as they they are tidally squeezed by 
the gravitational potential of the elliptical. Contrary to 
any model driven by major mergers, for which 1 Gyr old 
starbursts will be as stochastic as is the merging process, 
the rate of accretion should vary little with galaxy type, 
luminosity and perhaps even environment. 

Both scenarios ought to be tested quantitatively in the 
context of evolutionary models for the PN luminosities. 


stellar masses and galaxy spectra. The mass accretion 
rate can be obtained (Salvador-Sole et al. 1998) through 
the extended Press-Schechter formalism (e.g. Lacey and 
Cole 1993). Will there be sufficient blue stragglers in 
ellipticals to produce the universality of the PNLF? And, 
on the other hand, will the rate of accretion of small 
galaxies be large enough to produce sufficient bright PNe 
and not too large to make elliptical galaxies bluer than 
they are observed to be? 
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